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ABSTRACT 

Whether supernovae are major sources of dust in galaxies is a long-standing debate. We present 
infrared and submillimeter photometry and spectroscopy from the Herschel Space Observatory of the 
Crab Nebula between 51 and 670 /im as part of the Mass Loss from Evolved StarS program. We 
compare the emission detected with Herschel with multiwavelength data including millimeter, radio, 
mid-infrared and archive optical images. We carefully remove the synchrotron component using the 
Herschel and Planck fluxes measured in the same epoch. The contribution from line emission is 
removed using Herschel spectroscopy combined with Infrared Space Observatory archive data. Several 
forbidden lines of carbon, oxygen and nitrogen arc detected where multiple velocity components are 
resolved, deduced to be from the nitrogen-depleted, carbon-rich ejecta. No spectral lines are detected 
in the SPIRE wavebands; in the PACS bands, the line contribution is 5% and 10% at 70 and 100 /im 
and negligible at 160 /im. After subtracting the synchrotron and line emission, the remaining far- 
infrared continuum can be fit with two dust components. Assuming standard interstellar silicates, 
the mass of the cooler component is O.24^QQgM for T = 28.lig'2 K. Amorphous carbon grains 

require 0.11 ± 0.01 M of dust with T = 33.8l^'gK. A single temperature modified blackbody with 
0.14M Q and 0.08 M for silicate and carbon dust respectively, provides an adequate fit to the far- 
infrared region of the spectral energy distribution but is a poor fit at 24-500 /im. The Crab Nebula 
has condensed most of the relevant refractory elements into dust, suggesting the formation of dust in 
core-collapse supernova ejecta is efficient. 

Subject headings: dust, extinction-ISM: individual objects (Crab Nebula)-ISM: supernova remnants 
-submillimeter: ISM 



1. INTRODUCTION 

In galaxies, the major dust source has in the past 
been presumed to be low-intermediate mass stars during 
their asymptotic giant branch (AGB) phase, but when 
accounting for dust destruction timescales (e.g. Jones 
2001; Draine 2009) and the observed total dust masses, 
the required dust injection rate from stars can be an or- 
der of magnitude higher than observed (e.g. Matsuura et 
al. 2009; Gall et al. 2011; Dunne et al. 2011; Rowlands 
et al. 2012). An alternative source of dust is required to 
make up the dust budget (Pipino et al. 2011; Dunne et 
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al. 2011). This shortfall in the dust mass estimated from 
AGB stars is also observed in dusty high-redshift galax- 
ies where the timescales for dust production are close to, 
or shorter than, the lifetime of a typical low-mass AGB 
star (Morgan & Edmunds 2003; Dwek et al. 2007). 

Significant dust production in supernova (SN) ejecta 
would alleviate these dust budgetary problems. SNe have 
long been proposed as a source of dust (e.g. Dwek & 
Scalo 1980; Clayton et al. 2001) yet subsequent mid and 
far-infrared (FIR) observations have detected only small 
quantities of warm dust in young ejecta (Sugerman et al. 
2006; Meikle et al. 2011; Kotak et al. 2009; Andrews et 
al. 2011; Fabbri et al. 2011) and old remnants (Williams 
et al. 2006; Rho et al. 2008). These observed dust masses 
are orders of magnitude lower than required. 

In the era of the Herschel Space Observatory (Pilbratt 
et al. 2010), we are now piecing together the relative con- 
tribution of stellar sources to the dust budget in galax- 
ies, yet dust yields from the limited FIR studies of core- 
collapse remnants remain highly uncertain (Dunne et al. 
2003; 2009; Krause et al. 2004; Barlow et al. 2010; Mat- 
suura et al. 2011). Observations of historical Galactic 
remnants are important since these are (1) resolved, so 
that the different SN and interstellar/circumstellar trac- 
ers can be separated, and (2) young enough to ensure the 
thermal emission is not dominated by swept-up material. 
The Crab Nebula is one of only a few sources which sat- 
isfy these criteria and was chosen to be observed as part 
of the Herschel guaranteed time project Mass Loss from 
Evolved StarS (Groenewegen et al. 2011). This survey 
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includes Cas A (Barlow et al. 2010) and the Type la 
remnants Kepler and Tycho (Gomez et al. 2012; see also 
Morgan et al. 2003 and Gomez et al. 2009). Unlike these 
remnants, the Crab has negligible cirrus contamination 
along the line of sight and is an ideal source to minimize 
the effects of unrelated interstellar material. 

The Crab Nebula has been an object of interest for 
a number of years (see Hester 2008 for a comprehen- 
sive review). The remnant of an explosion in 1054 AD, 
the Crab is a pulsar wind nebula lying at a distance of 
2kpc (Trimble 1968). Its structure can be separated into 
two major components: the pulsar wind nebula (seen 
in X-rays and optical) with smooth synchrotron emis- 
sion (at near-IR and radio wavelengths), and a network 
of filaments (traced in the optical and IR) consisting of 
thermal ejecta. The low expansion velocity of the ejecta 
suggests the remnant is the result of a Type ITP ex- 
plosion (Mac Alpine & Satterfield 2008). This is further 
confirmed by abundance constraints which put the pro- 
genitor star at 9-12 M Q (Nomoto et al. 1982; Nomoto 
1985; MacAlpine & Satterfield 2008). 

Unusually amongst supernova remnants (SNRs), the 
material in the Crab is primarily photoionized by non- 
thermal radiation from the synchrotron nebula (e.g. 
Davidson & Fesen 1985; MacAlpine & Satterfield 2008). 
The latter authors found the main nebular gas compo- 
nent to be highly nitrogen-depleted and carbon-rich, al- 
though two other gas components with C/O < 1 were 
also found to be present. The implication of this is that 
both carbon-rich and oxygen- rich dust species could exist 
in the nebula. Using optical data, Woltjer & Veron-Cetty 
(1987) detected the presence of absorption attributable 
to dust at the position of a bright [O in] filament. Fesen 
& Blair (1990) obtained high angular resolution optical 
continuum images, revealing the presence of numerous 
"dark spots" across the synchrotron nebula coincident 
with bright emission cores seen in narrow-band [O i] , [C i] 
and [S n] images, consistent with dust existing in par- 
tially ionized or neutral clumps. 

Previous IR studies confirmed the presence of dust 
grains in the Crab as early as the 1980s with masses rang- 
ing from 0.005 to 0.03 Mq (e.g. Marsden et al. 1984). 
Green et al. (2004) used ISO and SCUBA to infer the 
presence of 0.02-0.07 M Q of dust depending on its com- 
position. A careful analysis of the line contribution to 
the mid-IR using Spitzer data out to 70 /zm suggested a 
dust mass of 0.001-0.01 Mq (Temim et al. 2006). Using 
spatially resolved MIR spectroscopy across the remnant, 
Temim et al. (2012) later revised this to a silicate grain 
mass of (2.4±f;|) x 10- 3 M Q . 

Previous studies either lacked long wavelength spec- 
troscopic information or adequate sampling of the FIR 
emission. In this paper, we present detailed FIR and cru- 
cially, submillimeter-through-to-millimeter observations 
of the Crab Nebula obtained with WISE, Spitzer, ISO, 
Herschel and Planck which allow us to accurately de- 
termine the synchrotron contribution and line emission 
to beyond 600 /im. We can then estimate the total dust 
mass within the ejecta for the first time. 



2. OBSERVATIONS AND DATA REDUCTION 
2.1. Herschel Photometric Imaging 
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Photometric Imaging 
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PACS 
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1342204441 
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PACS 


05 h 34 m 31 s .97, 22°00'52" 
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PACS 
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PACS 


05 h 34 m 31 s .97, 22°00'52" 
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ISO-3 


05 h 34 m 29 s .31, 22°00'37" 
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69301543 
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ISO-4 


05 h 34 m 34M9, 21°59'54" 
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69301611 


1630 s 


FTS 


05 h 34 m 29 s .47, 22°00'30" 


A 


1342204022 


3476 s 


IFU 


05 h 34 m 29 s .44, 22°00'32" 


.52 


1342217847 


2267 s 


IFU 


05 h 34 m 29 s .42, 22°00'47" 


.17 


1342217847 


1139 s 



Table lObservation Log for Herschel PACS and SPIRE Photom- 
etry and ISO LWS, Herschel PACS IFU, and SPIRE FTS Data Sets 
of the Crab Nebula. 




Figure 1. 3 X 3arcmin region showing the Herschel PACS image 
of the Crab Nebula at 70 /im with the locations of the ISO LWS 
apertures (#1-4 -small cyan circles), PACS IFU (white rectangle) 
and SPIRE FTS FoV (yellow dashed). The unvignetted FTS FoV 
is also shown (solid yellow circle) . Black crosses mark the locations 
of the Spitzer spectra from Temim et al. (2012). 

The Crab SNR was observed with the Herschel Pho- 
todetector Array Camera (PACS; Poglitsch et al. 2010) 
and Spectral and Photometric Imaging Receiver (SPIRE; 
Griffin et al. 2010) at 70, 100, 160, 250, 350 and 500 /im 
(a summary of the observations is listed in Tabled]). The 
PACS photometry data were obtained in "scan-map" 
mode with speed 20 arcsec s -1 including a pair of orthog- 
onal cross-scans over 22 arcsec x 22 arcsec. In order to 
obtain images at all three PACS wavelengths, we used 
both the 70+160 /im and 100+160 /jm channels, leading 
to the 160-/zm image having twice the exposure time of 
the other two channels. The SPIRE maps are "Large 
Map" mode with scan length of 30 arcmin over 32 arcmin 
x 32 arcmin at a speed of 30 arcsec s _1 ; a cross-scan is 
also taken, with a repetition factor of three. The data 
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were processed following the description given in Groe- 
newegen et al. (2011). 

The PACS photometric data were reduced with the 
Herschel Interactive Processing Environment (HIPE; Ott 
2010) applying all low-level reduction steps (including 
deglitching) to Level 1. The SCANAMORPHOS software 
(Roussel 2012) was then used to remove effects due to 
thermal drifts and uncorrelated 1// noise of the indi- 
vidual bolometers and create the Level 2 map. The full 
width at half maximum (FWHM) at 70, 100, and 160 /mi 
is 6, 8, and 12arcsec, respectively. The flux calibration 
uncertainty for PACS is less than 10% (Poglitsch et al. 
2010) and the expected color corrections are small com- 
pared to the calibration errors. We therefore adopt a 
10% calibration error. 

For SPIRE, the standard photometer pipeline (HIPE 
v.5.0) was used (Griffin et al. 2010) with an additional 
iterative baseline removal step (e.g. Bendo et al. 2010). 
The SPIRE maps were created with the standard naive 
mapper (e.g. Griffin et al. 2008). We multiply the 
350 /im data product by 1 .0067 to be in line with the most 
recent calibration pipeline (v7). The FWHM for pixel 
sizes of 6, 10, and 14arcsec is 18.1, 24.9, and 36.4arcsec 
at 250, 350 and 500 /im respectively. The SPIRE calibra- 
tion methods and accuracies are outlined by Swinyard et 
al. (2010) and are estimated to be 7%. The pipeline pro- 
duces monochromatic flux densities for point sources but 
at the longer wavelengths, color corrections become sig- 
nificant and we therefore use the correction factors listed 
in the SPIRE Observer's Manual (2011). 

2.2. Herschel PACS and SPIRE Spectroscopy 

A 51-210-/tm full-range spectral scan was obtained 
with the PACS Integral Field Unit (IFU) Spectrometer 
(Poglitsch et al. 2010). The IFU has 5x5 spaxels with 
each spaxel being 9.4 arcsec on a side. The coordinates 
(positioned on the brightest nebular filament - Figured]), 
are listed in TableQ] The "chop-nod" spectrometer mode 
was used, with the "off" position located 6 arcmin to the 
north and south of the "on" position for the red and blue 
range scans, respectively. The "chop-nodded" observa- 
tions contained one single nodding cycle and one single 
up-down scan in wavelength. The data were reduced 
to Level 2 products using the standard PACS chopped 
large range scan and spectral energy distribution (SED) 
pipeline in HIPE version 8.0.1 (Ott 2010) using calibra- 
tion file PACS_CAL_32_0. 

The SPIRE Fourier transform spectrometer (FTS) was 
used to obtain sparse-map spectra. Two bolometer ar- 
rays provided overlapping bands covering 32.0-51.5 cm -1 
(194-313 /tm for the SPIRE Short Wavelength Spectrom- 
eter Array-SSW) and 14.9-33.0 cm" 1 (303-671 /im for 
the SPIRE Long Wavelength Spectrometer Array-SLW). 
The SSW and SLW beamsizes are ~18 and ~37 arcsec, 
respectively. The two yellow circles in Figure [T] show 
the 2.2-arcmin diameter unvignetted and 3.2-arcmin di- 
ameter partially-vignetted field of view (FoV) for the 
SPIRE spectrometer pointing. A point-source calibra- 
tion was applied to the central detectors of each array 
using Uranus (Swinyard et al. 2010). As the Crab is es- 
sentially fully extended in the beam, the calibration is 
modified using the self-emission of the telescope as the 
primary source (as we know both its temperature and 
emissivity properties). This telescope Relative Spectral 



Response Function is photometrically calibrated against 
Uranus using knowledge of the instrument spatial re- 
sponse function as described in the SPIRE Observers 
Manual (2011). The resulting spectra have an absolute 
calibration accuracy of 5% compared to the Uranus flux 
model of G. Orton et al. (in preparation). 

2.3. ISO LWS Far-infrared Spectroscopy 

In addition to our Herschel spectroscopy, we have made 
use of archival 43-197-/xm spectra of the Crab obtained 
with the Long Wavelength Spectrometer (LWS; Clegg et 
al. 1996) aboard 750. Figure Q] shows the positions of 
the four pointings of the LWS superposed on the PACS 
70-/xm image of the Crab (cyan circles - Table [IJ. The 
44-110-/im parts of the spectra were previously published 
in Green et al. (2004). The data used here have been 
processed through the 75'0-LWS Highly Processed Data 
Products pipeline (Gry et al. 2003) and further modified 
by removing gain shifts between the individual detectors 
using the central 100-/im detector as the reference. The 
Crab is relatively faint and the continuum spectrum be- 
yond 150 /im suffers from poor dark current removal and 
possible non-linear response in the detectors due to ther- 
mal instabilities. This portion of the continuum is con- 
sidered untrustworthy although the line fluxes, in partic- 
ular the [Cn] 158-/tm line, are well calibrated. 

The global integrated fluxes at the Herschel wave- 
lengths are measured using an elliptical aperture with 
size 245 x 163" centered on the pulsar with the back- 
ground contribution estimated using apertures off the 
remnant. These fluxes are combined with IR and submm 
fluxes from the literature as described in detail in Ap- 
pendix [A] with photometric fluxes from 3.4 to 10,000 /tm 
listed in Table IA1I (see also Figures 0] and [5]) . These are 
combined with our own flux measurements at 3.4-22 /im 
with WISE and 3.6-70 /im Spitzer data (Appendix [A) . 

When comparing the IR fluxes measured in this work 
(marked as "A" in Table lATj) with the literature, we find 
that our PACS 70 /im flux for the Crab Nebula is sig- 
nificantly larger than the Spitzer measurement at the 
same wavelength (Temim et al. 2006). At first glance 
this appears to suggest a huge (40%) discrepancy be- 
tween PACS and MIPS measurements (e.g., Aniano et 
al. 2012). However, we find this discrepancy disappears 
if we use the most recent calibration factors from Gordon 
et al. (2007). Re-reducing the 70 /im map using the most 
recent DAT instrument team pipeline produces a 70 /im 
flux of 210 Jy after the relevant color corrections are ap- 
plied, this is in excellent agreement with our PACS mea- 
surement. Note that claimed differences between PACS 
and Spitzer MIPS can be resolved if instrumental ef- 
fects are carefully considered; the PACS detectors are 
extremely stable with virtually no non-linearities com- 
pared to the Ge:Ga photoconductors employed in MIPS 
and indeed ISO. For these reasons, PACS has achieved 
an absolute calibration accuracy of 5% for point sources 
with < 10% quoted here (compared to 20% for MIPS at 
70 /im). 

Since the radio synchrotron flux decreases with time 
by -0.167 %yr _1 (Aller & Reynolds 1985), radio fluxes 
measured at different epochs (Table lATj) need to be cor- 
rected to the 2010 epoch to allow for comparison with the 
WISE, Herschel, and Planck data. Note that this "fad- 
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Figure 2. Multiwavelength montage of the Crab SNR centered on the pulsar (R.A. = 262°. 671, decl. = 22°. 0145, J2000.0) with diameter 
4'. 2 X 4'. 2. Top: (from left to right): three color image in Ho (red), [Oiii] (green) and Bessel B (blue) using images obtained from the 
Faulkcs Telescope North (Appendix [Aj . Spitzer IRAC 4.5 (im, MIPS 24 /im and VLA archive data at 1.4 GHz (taken in 1996). Middle: 
Herschel PACS three color image and individually the 70, 100 and 160 fim images. Bottom: Herschel SPIRE three color image and 
individually, 250, 350 and 500 (im images. Beam sizes are indicated with the black circle in the lower left corner. 



ing rate" assumes that the non-thermal fluxes decline at 
the same rate at all wavelengths. 

3. DISENTANGLING THE CONTRIBUTIONS FROM 
DIFFERENT COMPONENTS 

In Figure [H we present a multiwavelength view cen- 
tered on the Crab, comparing the Herschel PACS and 
SPIRE images with synchrotron emission seen at near- 
IR (with Spitzer IRAC) and at radio wavelengths with 
the Very Large Array (VLA), and the ionized gas ob- 
served at optical wavelengths. The smooth non-thermal 
synchrotron emission seen at 4.5 /im appears to be con- 
fined within the filamentary structures seen in optical 
and IR wavebands which originate from ionized gas and 
dust emission. Previous works have shown that the warm 
dust component (seen in emission at 24 /im) traces the 
densest gas in the cores of filaments in low ionization 
states (e.g., [O i] - Blair et al. 1997; Loll et al 2007) and 
this is confirmed by the filamentary emission seen in the 
MIPS and Herschel PACS images. The distribution of 



the emission in the Herschel SPIRE 250 /im image is sim- 
ilar to the shorter IR wavelengths. At the longer SPIRE 
wavelengths (Figure [2]), we start to see a strong resem- 
blance with the radio emission at 1.4 GHz, which traces 
both the smooth synchrotron seen at the shorter 4.5 /im 
band, but also some filamentary emission arising from 
free- free emission (e.g. Temim et al. 2006). While the 
free-free emission only makes a negligible contribution to 
the FIR continuum for the Crab (see Appendix IB)) . the 
synchrotron component and line emission are important 
at these wavelengths and need to be removed before we 
can investigate if there is residual emission from dust. 

3.1. Line Emission 

Temim et al. (2012) present a comprehensive analy- 
sis of the MIR spectral lines across the Crab Nebula, 
with a number of forbidden lines identified up to 36 /im. 
They estimate that the contribution of line emission to 
the 24 /tm broadband flux is 27%, 54%, and 48% mea- 
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Figure 3. Velocity profiles of the emission lines detected in the Crab Nebula using the coadded PACS spectra (red) and the LWS spectrum 
(black) taken at the same location (#3). The text in the top right corner lists the approximate central wavelength of the lines which 
at 51.8145 /xm and 88.3564 fim for [O III] (NIST Atomic Spectra Database); 63.18371/mi and 145.5255 (im for [O i] (Zink et al. 1991); 
157.74095 /im for [C II] (Cooksy et al. 1986); and 121.8976 /im and 205.1783 fim for [N n] (Brown et al. 1994). 



sured at different locations across the remnant (shown 
in Figure Q}, suggesting on average, that line emission 
contributes 43% ± 6% of the flux (Tabled]). (The error 
quoted in this average is simply the range of values ob- 
tained on the dense filaments.) Note that this estimate 
of the line contribution was obtained after Temim et al. 
subtracted synchrotron emission from the Spitzer map. 



A 


Line Contribution in 


Band (%) 


(fan) 








Average 


24 




54 48 


27 


43 ± 6 a 




#1 


#2 #3 


#4 




70 


6.6 


4.2 5.4 


4.4 


4.90 ± 0.05 b 


100 


12.9 


5.6 9.6 


6.7 


8.7 ± 0.3 b 



Table 2The contribution to the total integrated flux at 24, 70 
and 100 fim at different locations across the remnant (Figure [TJ- a 
- This represents the contribution of line emission to the flux after 
synchrotron subtraction has been carried out (Temim et al. 2012). 
b - The contribution at 70 and 100 fim from line emission estimated 
using ISO LWS and Herschel PACS. 

In order to determine the contribution of line emission 
to the broad band infrared fluxes of the Crab Nebula, 
Herschel PACS and SPIRE, and ISO LWS data were an- 
alyzed. Figure [3] shows the velocity profiles of the seven 



emission lines detected in the co-added Herschel PACS 
spectra and in the LWS spectrum taken at the same lo- 
cation. Although the LWS line profiles are significantly 
broadened compared to the instrumental resolution, it 
requires the higher PACS spectral resolutions to fully re- 
solve the individual blueshifted and redshifted velocity 
components from the ejecta. The separations between 
the blue and red component emission peaks vary between 
1290 and 1750kms _1 , depending on the species (the pro- 
files will be discussed further in a separate paper). 

The seven emission lines and their gas properties are 
discussed fully in Appendix [Cl with flux es m easured rel- 
ative to the [Om] line provided in Table [CT1 These lines 
suggest the Crab Nebula ejecta is nitrogen-depleted and 
carbon-rich (unlike Cas A; Rho et al. 2008) implying that 
the dust is likely to be carbon-rich. However, there is also 
evidence that some regions in the Nebula have solar-like 
CNO abundances which may allow silicate- type grains to 
also form (see Appendix [Cl for more details). 

The four LWS spectra were used to estimate the contri- 
butions from the emission lines to the measured PACS 70 
and 100 /im broad-band fluxes. To do this, the PACS fil- 
ter spectral response functions were convolved with and 
integrated across (1) the observed LWS spectra, includ- 
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Figure 4. IR SED of the Crab Nebula including the integrated 
fluxes from Herschel, the new Spitzer calibration and WISE fluxes 
(black points). The errors inclu de th e photometric errors (domi- 
nated by calibration errors, Table EU, the error in the synchrotron 
extrapolation and the error in the line emission contribution (Ta- 
ble [2$ added in quadrature. The ISO #3 spectra scaled to the 
100 (im flux is overplotted in gray. The average (dereddened - 
see Indebetouw et al. 2005) Spitzer spectroscopy (from the two 
brightest filaments in Figure [TJ is overplotted for comparison. The 
synchrotron power law is shown in red. Two temperature compo- 
nents of amorphous carbon dust at 63 and 34 K required to fit the 
SED are plotted (blue dashed) with the sum of these plotted in 
solid blue. A single-temperature amorphous carbon fit (40 K) to 
the FIR is also shown (green dot-dashed). The sum of the two- 
component dust and synchrotron contributions is shown by the 
solid black curve. 

ing emission lines; and (2) the LWS spectra after ex- 
cising emission lines and interpolating the adjacent con- 
tinua across the positions of the lines. The mean ra- 
tio of the case (2) to case (1) in-band fluxes was found 
to be 0.951±0.012 for the 70 /im filter and 0.913±0.028 
for the 100 /im filter (Table [5]). We have applied these 
average line correction factors to the measured 70- and 
100-^m fluxes to obtain the continuum-only broad band 
fluxes. The 146 and 158 /im emission lines (Table IClf) 
contribute a negligible amount to the broad-band PACS 
160 /im flux measurement. No emission lines were de- 
tected in any of the spectra from the SPIRE FTS detec- 
tors, the contribution of line emission to the broad-band 
SPIRE photometry is therefore negligible. 

3.2. Synchrotron Emission 

In this section, we now estimate the contribution from 
the non-thermal synchrotron component. To do this, we 
fit a power law to the Spitzer IRAC, WISE and Planck 
data (large circles in Figure [5]). 

Including the Planck fluxes in the SED allows us 
to fully constrain the synchrotron power law between 
3.4 /im and 1cm for the first time, using observations 
taken at the same epoch. A least-squares fit to the 
Spitzer- WISE- Herschel- Planck data set produces a power 
law with frequency dependence i/~ 0A17 and amplitude 
1489 Jy at 1 GHz (Figured]). The error in the syn- 
chrotron spectral index (a) from the line of best fit is 
±0.006. Extrapolating the expected fluxes at the IR- 
submm wavebands due to synchrotron emission using the 
integrated fluxes for the remnant with the above power 
law, we estimate the synchrotron flux in each waveband 
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Figure 5. SED of the Crab Nebula from the IR-radio including 
Herschel (red points) and Spitzer (black), WISE (purple) photom- 
etry and Planck fluxes (Planck Collaboration 2011 , blue points). 
Previous fluxes from the literature (see Table IA1I and references 
therein) are shown with gray diamonds. The synchrotron law fit- 
ted to the 3.6 — 10 4 /im data points is the dashed black line. Dot- 
dashed lines are the fitted components from thermal emission by 
amorphous carbon grains (see Figure |4jl. The solid black line is 
the total flux obtained from summing the synchrotron and the two 
dust components (with the residual shown below). Note the total 
integrated fluxes plotted here also includ e a c ontribution from line 
emission at 24, 70, and 100 fim (Section 13.111 . which has not been 
added to the total black SED curve plotted here. 

(Table EE]). 

The synchrotron power law in the FIR regime derived 
in this work is steeper than the slope determined in the 
review by Macfas-Perez et al. (2010), where the extrapo- 
lation to the submm from the longest radio wavelengths 
suggested one synchrotron component with a = —0.3. 
We argue here that the exquisite coverage between the 
IR and radio regime given by Herschel SPIRE and Planck 
suggests that the synchrotron emission is described by (at 
least) two power laws, with wavelengths beyond 10 4 /im 
(30 GHz), following the flatter relationship i/ -0 3 (as also 
seen in Green et al. 2004; Arendt et al. 2011). The slope 
then breaks, becoming the steeper 2 law we find 
here. In fitting this synchrotron power law, we do not 
incorporate the previous literature data obtained in the 
submillimeter-millimeter regime (shown by the faint gray 
diamonds in Figure [5]) since these fluxes often have large 
calibration errors, with small FoVs. One should be par- 
ticularly cautious when using the ISOPHOT data of the 
Crab (Green et al. 2004) since this was obtained in P32 
chopping mode, which can suffer significantly from tran- 
sient effects with unreliable calibration (U. Klaas, pri- 
vate communication). This dataset was never released 
as a scientifically validated measurement and therefore 
the quoted calibration accuracy for ISOPHOT (30%) is 
not applicable for this dataset. The previous literature 
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measurements were also taken at different epochs and 
corrected using the average "fading" rate (Appendix [A"|l . 
With WISE, Herschel, and Planck, not only are the pho- 
tometric errors less than 10%, the images also have a 
large background area to sample and the data were taken 
at the same epoch, therefore not relying on the applica- 
tion of a correction factor. 

It is possible that the synchrotron spectrum could fur- 
ther break into different components in the IR regime (as 
suggested by Arendt et al. 2011) which would introduce 
further errors on the amount of synchrotron estimated at 
each wavelength. However, we see no evidence for a sharp 
break either via the presence of excess flux or in differ- 
ences in the spectral index maps created from the IRAC 
bands versus maps created from the VhA-Herschel im- 
ages. Unfortunately, given the low angular resolution in 
the FIR/submm, any local variations in the spectral in- 
dex would not be seen in the method we have used in this 
work. The present data are therefore insufficient to sep- 
arate out any small-scale variations in the synchrotron 
slope which may account for some of the residual emis- 
sion. The shape of the SED at wavelengths beyond 70 /im 
and the Planck coverage rules out a significant break in 
the FIR/submm regime at least. However, we note that 
the flux attributed to synchrotron could be underesti- 
mated at 24 /im if there is a break to a steeper power law 
at wavelengths less than 70 /im. 

In order to spatially determine the distribution of syn- 
chrotron and excess thermal emission from dust, we fol- 
low Temim et al. (2006) who used a combination of 
the Spitzer IRAC and radio images to obtain a spec- 
tral index map. We repeat this procedure using the 
extinction-corrected 4.5 /im IRAC image and the 500 /im 
Herschel SPIRE map, both of which are completely dom- 
inated by synchrotron emission. After aligning and con- 
volving the IR-submm images, we re-grid them to the 
same pixel size and deconvolve our spectral index im- 
age with the appropriate beam. We use the deconvolved 
spectral index map to subtract the extrapolated syn- 
chrotron flux expected at each IR-submm waveband us- 
ing the spectral index for that pixel. In Figure|6]we show 
the extrapolated synchrotron emission expected at 24 /im 
using the spectral index map made in this way. Note that 
the subtraction of the synchrotron component on a pixel- 
by-pixel basis removes the smooth emission at 24 /im seen 
in Figure [2j leaving only the excess filamentary emission 
originating from the warm dust component. 

3.3. Thermal Emission from Dust 

Using the global SED of the Crab (Figures 0] and [5]), 
after removing the contribution from line emission and 
the well-constrained synchrotron component, we find the 
excess thermal emission observed in the FIR can be de- 
scribed by the sum of two modified blackbodies arising 
from a warm (T w ) and a cool (T c ) component. Although 
we would expect a more complex temperature distribu- 
tion for dust in the remnant, a two-temperature compo- 
nent fit is adequate for a first-order approach in modeling 
the SED. We can then fit the data with optical constants 
appropriate for silicate or carbon grains, with the dust 
mass estimated using Equation ([I}. S v is the flux den- 
sity, D is the distance, B{y,T) is the Planck function 
and k„ is the dust mass absorption coefficient calculated 
from the dust emissivity Q u (for silicates - Weingartner 



& Draine 2001), and the grain density p (Laor & Draine 
1993): 



The total dust mass from the FIR model (see Tabled 
for astronomical silicates is dominated by the cool com- 
ponent which requires a best-fit temperature of T c — 
28.1+!;! K and mass M d = O.24±°;^M . The warm 
component (Figure 2]) arises from (8.3±|;|) x 10~ 3 M Q of 
dust at a temperature of 55.6+2 sK. 

As the gas in the filaments is carbon-rich (Ap- 
pendix fUJi the dust may be composed of amorphous car- 
bon grains. In this case (with Q v taken from Zubko et al. 
1996 (their 'BE' model) and p from Rouleau & Martin 
1991), the dust mass arising from the cool component re- 
duces to 0.11 ± 0.01 M Q with T c = 33.8+^K. The warm 
component requires only Md = (6.0+ 2 ' 4 ) x 10~ 3 M Q at 
T w = 63.4+J^K. 

We also attempted to fit iron grains to the SED (fol- 
lowing Matsuura et al. 2011) yet it is difficult to fit the 
MIR part of the thermal emission. To explain the FIR 
emission (i.e. beyond 70 /im) with iron grains, we require 
~ 0.3 Mq of dust at temperatures 34 and 69 K with ra- 
dius 0.1 /im (where Q v is taken from Semenov et al. 2003 
and p from Nozawa et al. 2006). However, the mass of 
iron grains is highly dependent on the grain size, for ex- 
ample grains with radius < 0.005 /im would require more 
than 70 M Q of dust. 

We attempted to fit the residual emission from 24 
to 350 /im with a single-temperature modified black- 
body, producing the parameters Md ~ 0.14 Mq with 
T ~ 34 K for astronomical silicates and Md ~ 0.08 Mq 
with T ~ 40 K for amorphous carbon (the green dot- 
dashed curve in Figure [J). However, the single compo- 
nent always severely underestimates the flux at 24 /im. 
For this to be explained via contamination of line emis- 
sion in this waveband, we would require 96% of the 
broad-band flux at 24 /im to be due to line contribution; 
this is clearly not supported by the careful analysis of IR 
spectra across the remnant in Temim et al. (2012). Given 
that the two-component model adequately fits the entire 
IR-submm SED and that the chi-squared statistic favors 
a two-component model fit (Table [3]), we suggest this is 
the most valid model for the data. Ultimately, the differ- 
ence in the final mass is small whether one or two dust 
populations exist if the grains are carbon-rich though the 
silicate dust mass estimated from the single-temperature 
fit is reduced by approximately half. 

Note that modeling the SED with dust grains with 
a continuous temperature distribution could produce a 
lower dust mass compared to that estimated using the 
two-temperature approach. The dust masses quoted 
in Table [3] should then be regarded as an upper limit 
on the mass of dust in the remnant. However the 
largest uncertainty in a quoted dust mass arises from 
the choice of optical constant. This becomes more 
important when choosing models within the "umbrella" 
term amorphous carbon, which encompasses a number 
of different classes of materials such as soot, glassy 
carbon (Jageret al. 1998), and the ACAR/ACH2/BE 
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Two-component Model 
Warm Component Cool Component 
T W (K) M d (xlO- 3 M ) T c (K) M d (M ) 


X 2 


Si 


55.61™ 




28.i±i;i 


°- 24 -o.i 


0.05 


c 


63.4±f;* 


e.o±J:i 


33.8± 2 ;| 


0.11 ±0.01 


0.23 




T d (K) 


One-component Model 

M d (M ) x 2 






Si 
C 


31 
40 


0.14 
0.08 


3 x 10~ 5 
7 x 10~ 3 







Table 3Summary of the best-fit dust masses and temperatures 
for silicate and amorphous carbon grains using the two (top) and 
one (bottom) component dust models displayed in Figures [4] and 
\E\ Note: the reduced x 2 statistic is also included. 



models of amorphous grains from Zubko et al. (1996). 
Although the optical constants vary depending on which 
of these grain types is chosen, at FIR wavelengths Q v 
varies by only a factor of a few for these different types 
(Figure 10 in Jager et al. 1998; Hanner 1988). One ex- 
ception to this arises if the grains are pyrolyzed cellulose 
where Q„ is different by a factor of 10 (Jager et al. 1998). 

The distribution of the warm dust component is shown 
in Figure [6] using the synchrotron- and line emission- 
subtracted 24 /im image. We can now use this warm 
dust map to extrapolate the expected emission from the 
warm dust component in the Herschel bands and reveal 
the distribution of the cool dust component. The spatial 
distribution of the synchrotron emission at 160 /im (using 
the spectral index map) is shown in Figure [6] along with 
the extrapolated warm and cool dust components at this 
wavelength. The morphology of this cool dust is similar 
to the warm dust emission at 24 //m and is distributed in 
the dense filamentary structures as expected. 

4. DUST IN THE CRAB NEBULA 

Using Spitzer data out to 70 /mi, Temim et al. (2012) 
derived a silicate grain mass of 2.4 x 10~ 3 M© (T ~ 55 K) 
and 3.2 x 10~ 3 M© (T ~ 60 K) for carbon grains. Their 
estimates are of similar order to the mass of grains 
estimated from the hot component listed in Table [3l 
However, the total dust mass derived here is at least 
an order of magnitude higher than Temim et al. with 
M d ~ 0.24 M @ (T c ~ 28 K) or 0.11 M© (T c ~ 34 K) 
for silicates and carbon grains. The difference in mass 
between these two studies is not due to the optical con- 
stants assumed, since Temim et al. use the ACAR model 
from Zubko et al. (1996) which has similar Q v to the 
BE model used here at FIR wavelengths. The differ- 
ence is due to the longer wavelength coverage of Her- 
schel PACS and SPIRE, the matched-epoch observations 
and careful subtraction of synchrotron emission in the 
FIR-submillimeter regime. 

To account for their dust temperature of 50-60 K, 
Temim et al. (2012) modelled the dust heating (via non- 
thermal radiation) and cooling rates (through IR emis- 
sion) in the Nebula. They showed that the IR emis- 
sion (A < 70 ^m) originates from small grains with radii 
< 0.015 /im. Their results imply that the SN grains in the 
Crab are small, and therefore easier to destroy via sput- 
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Figure 6. Top: the distribution of synchrotron and thermal emis- 
sion from the warm dust component at 24 /im. Bottom: the Her- 
schel SPIRE 160 (im emission separated into the synchrotron com- 
ponent (left), the warm dust component (as traced by the emission 
at 24 fim) (middle) and the newly identified cool dust component 
(right). The color palette used in the online version of this figure 
is the cubehelix scheme in which color monotonically increases in 
terms of perceived brightness (Green 2011). 

tering; the authors also point out that this is somewhat 
at odds with theoretical models of SN grain formation 
(e.g. Kozasa et al. 2009). Given the newly detected cool 
dust component in this work, the long- wavelength radia- 
tion originates from larger grains with radius > 0.06 /im. 
Such grains are predicted by the Kozasa et al. model and 
would suggest this new cool component of dust would be 
harder to destroy. 

Since the dust is spatially coincident with ionized 
ejecta material (Figure [3]), it is plausible that the ther- 
mal emission arises from newly formed grains. Indeed, we 
can rule out a swept-up interstellar origin using a sim- 
ple argument (see, for example, the Tycho and Kepler 
remnants-Gomez et al. 2012). The volume swept up by 
the Crab is ~ 15 pc 3 , sweeping up a total gas mass of 
approximately one-tenth of a solar mass (Trimble 1970; 
Davidson & Fesen 1985) for typical interstellar densi- 
ties. Applying a standard gas-to-dust ratio (Devereux & 
Young 1990), the swept-up dust mass would therefore be 
< 1O _2 M . At the location of the Crab (180 pc above 
the Galactic plane), the interstellar density is thought to 
be much less than the canonical value, indeed the sur- 
rounding medium appears devoid of material (Davidson 
& Fesen 1985), placing even more stringent constraints 
on the possibility that the dust here is swept up material. 

Are the dust masses estimated here sensible given the 
expected heavy element abundances from the SN? The 
total amount of heavy elements (and therefore the mass 
available to form dust) in the SN ejecta from a 9-12 M© 
progenitor star varies from 0.2 to 0.5 M© depending on 
the model used and the mass of the progenitor (Maeder 
1992; Woosley & Weaver 1995; Limongi & Chieffi 2003; 
Nomoto et al. 2006). The limit on the abundances ex- 
pected in the ejecta from these theoretical models pro- 
vides the constraint for the maximum possible dust mass, 
these are 0.09 M© for carbon, 0.03 M© for MgO, and 
0.04 M© for Si02- For silicate-rich grains, we can in- 
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crease the available mass to 0.19 Mq if we allow iron to 
form grains, for example including grain compositions 
such as FeSi03. 

From the elemental abundances constraint, both the 
silicate and carbon dust masses estimated from the SED 
fitting are well within the maximum possible dust masses 
allowed for these compositions. Whether the grains are 
silicate- or carbon-rich, the observed dust masses suggest 
efficient condensation of metals in the filaments. 

5. CONCLUSIONS AND DISCUSSION 

The combination of Spitzer, WISE, ISO, Herschel, 
Planck photometry and spectroscopy at longer wave- 
lengths than probed before, reveals a previously unknown 
cool dust component in the Crab Nebula located along 
the ionized filamentary structures. 

To reveal this new dust component, we carefully re- 
moved the synchrotron component using WISE, Her- 
schel and Planck fluxes for the remnant measured in 
the same epoch. We find a steeper V - QA17 power-law 
variation with frequency for the non-thermal component 
which describes the emission from synchrotron at wave- 
lengths of 3.4-10,000 urn. 

The contribution from line emission is then removed 
using Spitzer IR spectroscopy (Temim et al. 2012) and 
our analysis of Herschel spectroscopy combined with ISO 
LWS archive data. FIR spectroscopy yields high O/N 
ratios and although appears to be carbon-rich, also has 
a component with solar-like CNO abundances. It seems 
likely that carbon and silicate grains could be located in 
the ejecta. 

The mass of dust estimated using the two-temperature 
approach to fitting the SED ranges from 0.1 to 0.2 Mq 
depending on whether the grains are composed of amor- 
phous carbon or astronomical silicates respectively. The 
warm and cool dust components are distributed in the 
filaments coinciding spatially with Doppler-shifted ejecta 
material (traced via spectroscopy). This indicates the 
dust is spatially coincident with the ejecta material and 
not swept up circumstellar/interstellar material. 

Comparing with the expected elemental abundances in 
the ejecta, this work suggests the condensation of heavy 
elements into dust grains is efficient, and that the fila- 
ments may provide a viable environment to protect the 
dust from shocks. The dust mass estimated for the Crab 
in this work is similar to the cool (unambiguously) as- 
sociated SN dust mass observed in Cas A, and at the 
lower limit of the mass estimated in SN1987A using re- 
cent Herschel observations (Matsuura et al. 2011). To- 
dini & Ferrara (2001) and Kozasa et al. (2009) predict 
between 0.1 and 0.3 Mq of dust should form in the ejecta 
from the (Type-IIP) explosions of progenitor stars with 
initial mass < 15 Mq, in agreement with the dust masses 
derived here. 

It is unclear how much of the newly formed ejecta dust 
will survive. Grains will be destroyed via thermal sput- 
tering in the reverse shock due to collisions with elec- 
trons and/or ions, yet unlike the Cas A and SN1987A 
remnants, the Crab does not have a visible reverse shock 
today (Hester 2008). The current environment in which 
the dust particles find themselves in appears relatively 
benign, as shown by the large amounts of molecular hy- 
drogen comfortably surviving within the filaments (Gra- 
ham et al. 1990; Loh et al. 2011). The dust particles in 



the Crab Nebula appear well set to survive their journey 
into the interstellar medium and contribute to the inter- 
stellar dust budget. Future ALMA observations of this 
source and other remnants are crucial to disentangle syn- 
chrotron and dust emission on smaller scales. This will 
be particularly important in comparing the time evolu- 
tion of dust forming and being destroyed in remnants at 
different stages. 

APPENDIX 
A. ANCILLARY DATA 

We used 3.6-70 /xm data from Spitzer IRAC and MIPS 
(PI. R. Gerhz; Temim et al. 2006). The fluxes were cor- 
rected for extended emission and color correction. Cal- 
ibration uncertainties were assumed to be 5% for IRAC 
and 10% at 24 and 70 fim and we applied the most up- 
to-date calibration factors from Gordon et al. (2007). 

We obtained single exposure (level lb) images of the 
WISE (Wright et al. 2010) all-sky release through the 
NASA/IPAC Infrared Science Archive at 3.4, 4.6 and 
22 fim. Montage was used to process and co-add the 
single exposure images. Median filtering of the single 
exposure frames was used to make basic cosmetic correc- 
tions. Calibration factors were applied according to the 
WISE Explanatory Supplement with photometric errors 
assumed to be 5% (see e.g Mainzer et al. 2012). 

For comparison purposes, we use optical images of the 
Crab (Figure[2]) using the 2.0 m Faulkes Telescope North 
on in Ha (200 s) and [O in] narrowband filters (240 s) and 
the broad-band Bessel blue filter (200 s). 

Photometric fluxes measured in previous works were 
added to the Herschel and Spitzer datasets presented 
here, including fluxes from Planck (Planck Collaboration 
2011), ARCHEOPS (Maci'as-Perez et al. 2010), theKuiper 
Airborne Observatory (Wright et al. 1979), IRAS (Strom 
& Greidanus 1992), ISOPHOT and SCUBA (Green et al. 
2004). The fluxes listed in Table ED that were not mea- 
sured in this work, were taken from the compilation of 
literature fluxes in Maci'as-Perez et al. (2010), Arendt 
et al. (2011), or from the Planck archive. In this work, 
we adopt calibration errors of 15% for IRAS. For the 
ISOPHOT data, we can only assume a calibration er- 
ror of 30% which is appropriate for scientifically valid 
datasets (from mode P22, for example). However, the 
Crab data have not been scientifically validated so this 
is a lower limit on the flux error. 

B. CONTRIBUTION FROM FREE-FREE EMISSION 

Using a simple argument we can demonstrate that 
free-free radiation makes a negligible contribution to the 
200 -/im continuum flux of the Crab Nebula. 

An upper limit can be derived starting from observed 
F(R(3) values (taken here to be 1.78 x 10~ n ergcm~ 2 s _1 
from MacAlpine & Uomoto 1991). They also measure a 
global i^(5876A)/F(H/3) ~ 0.84, while Davidson (1987) 
measured F(4686A)/F(H/3) - 0.74 and Smith (2003) 
measured ~ 0.72 giving a mean ratio of 0.73. De- 
reddening these ratios with E(B — V) = 0.47, adopting 
a mean temperature (T e ~ 9000 K) and using standard 
recombination emissivities, gives n(He + )/n(H + ) = 0.44 
and n(He 2+ )/n(H+) = 0.65. 
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Table A 1 Integrated fluxes for the Crab Nebula. 

Notes: Also included are the synchrotron fluxes derived from 
^-0.417 

References: A - this work; a - Temim et al. (2006); b - Green 
et al. (2004); c - Strom & Greidanus (1992); d - Wright et al. 
(1979); e - Arendt et al. (2011) /- Planck Collaboration 2011; g - 
Macfas-Perez et al (2010); h - Chini et al. (1984); i - Werner et al. 
(1977). * - These errors are the uncertainty in flux quoted in the 
Planck catalog which does not include calibration errors (< 7%). 
CH - University of Chicago photometer. KAO - Kuiper Airborne 
Observatory. 

Using the above numbers with Equation 6 from Milne 
& Aller (1975), substituting in F„(5GHz)/J(H/3), where 
F„(5GHz) is the 6cm free-free flux and /(H/3) = 8.6 x 
10 -11 ergcm~ 2 s _1 is the de-reddened HP flux, predicts 
a value at 5 GHz of F v = 0.88 Jy. Allowing for the i/~ - 09 
Gaunt factor frequency dependence of free-free emission, 
the extrapolated flux at 200 /im due to the free-free com- 
ponent is F U (E) = 0.53 Jy and 0.44 Jy at 24 /im. 

C. LINE EMISSION, ELECTRON DENSITIES, AND IONIC 
ABUNDANCES 

Table [Cll lists the integrated fluxes in the 88 /mi [O in] 
line measured in each of the four LWS spectra, as well 
as its flux in the co-added LWS spectrum and in the co- 
added spectra from the 25 PACS IFU spaxels. The table 
also lists the relative fluxes for the other detected lines 



(Figure GO), on a scale where F(88 /im) = 100.0. The 
[Nin] 57 ^m line is weakly detected in the spectra. The 
[O in] 52 yitm line shows a well-resolved double-peaked 
profile in the co-added PACS spectrum; however as it lies 
at the extreme short-wavelength end of the PACS wave- 
length coverage, where the responsivity is falling steeply, 
its relative flux calibration is uncertain and so we rely on 
the LWS measurements. 

Since the observed [O i] and [C n] lines can arise from 
mainly, or partly, neutral regions, we confine our anal- 
ysis here to lines that originate from mainly ionized re- 
gions of the nebula. The flux ratio of the [O III] 52 and 
88 /im lines is an electron density diagnostic (see Fig- 
ure 3 of Liu et al. 2001). We have used the same 2+ 
atomic datasets as Liu et al. to derive electron densi- 
ties of 135-485 cm" 3 from the LWS 52/88 line flux ratios 
listed in Table ICll These values are somewhat smaller 
than the values of 830-1230 cm -3 derived from the ra- 
tios of the [Sll] 6716 A and 6731 A line fluxes measured 
at several positions in the Nebula by MacAlpine et al. 
(1996), and from the ratios of the [S in] 18.7 and 33- 
5 /im line fluxes measured at several positions by Temim 
et al. (2012). HST imagery presented by Sankrit et al. 
(1998) has demonstrated that the optical [Om] emission 
from the Crab originates from diffuse sheaths around the 
filamentary cores that are bright in emission lines from 
lower ionization species, consistent with the lower [O in] 
densities that we find here. 

In Table [Cll we present flux ratios and the resulting de- 
rived ion ratios estimated using our tabulated n e (0 III) 
values and the atomic data sets adopted by Liu et al. 
(2001). Liu et al. (2001) have shown that, because of 
their similar critical densities, the ratio of the fluxes in 
the [O III] (52+88) /im lines to that in the [N in] 57 /tin 
line yields 2+ /N 2+ ion abundance ratios that are in- 
sensitive to the adopted electron density. In addition, 
because of their similar ionization potentials, 2+ /N 2+ 
is a good approximation to O/N. F(52 + 88)/.F(57) ra- 
tios measured from the LWS spectra, together with the 
Q 2 + /N 2+ ion ratios derived from them are listed in Ta- 
ble [Cll The flux ratio F(122)/F(57) and the resulting 
derived N + /N 2+ N + /N 2+ ratios show singly ionized ni- 
trogen to be the dominant nitrogen ion in the Crab. The 
2 +/N 2+ (O/N) ratios from the individual LWS spec- 
tra straddle a range of 10-34 by number. These values 
can be compared with the elemental abundances esti- 
mated for three 'Domains' by MacAlpine & Satterfield 
(2008). If we assume the solar abundances of Asplund et 
al. (2009) then using their Table 2, their mass fractions 
correspond to O/N number ratios of 7.2 (solar), 21 and 
260 for Domains 1, 2, and 3, respectively. The individual 
O/N values derived here from LWS spectra #1-3 are con- 
sistent with Macalpinc & Sattcrficld's Domain 2. They 
found this to be the most prevalent material, with nitro- 
gen depleted by a factor of three and carbon enhanced 
by a factor of six, relative to solar, which they noted to 
be consistent with a precursor star Mi > 9.5 Mq. 

We thank the referee Tea Temim for her construc- 
tive and helpful referee report. This research made use 
of APLpy, an open-source plotting package for Python 
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240 ± 30 


485 ± 30 


220 ± 20 


350 ± 30 


135 ± 20 




F(52 + 88)/F(57) 


18.5 ± 1 


18.9 ± 0.7 


42.8 ± 3 


12.8 ± 1 






2 +/N 2 + 


15.0 ± 0.9 


14.4 ± 0.5 


35 ± 3 


17 ± 1 






F(122)/F(57) 


0.29 ± 0.05 


0.54 ± 0.1 


0.51 ± 0.1 


0.16 ± 0.04 






N+/N 2 + 


2.3 ± 0.4 


5.3 ± 1 


3.9 ± 1 


1.4 ± 0.4 







Table CI Far-infrared line fluxes for the Crab Nebula including electron densities and relative ion abundances of the gas. The relative 
fluxes are given on a scale where _F(88 fim) = 100. rma The PACS relative flux calibration is uncertain at 52 um. 
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